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The Grover quantum search algorithm is generalized to deal with an arbitrary mixed initial 
state. The probability to measure a marked state as a function of time is calculated, and found 
to depend strongly on the specific initial state. The form of the function, though, remains as it 
is in the case of initial pure state. We study the role of the von Neumann entropy of the initial 
state, and show that the entropy cannot be a measure for the usefulness of the algorithm. We give 
few examples and show that for some extremely mixed initial states (carrying high entropy), the 
generalized Grover algorithm is considerably faster than any classical algorithm. 



m 
O 
o 

C 



> : 

m . 
oo ■ 

vo : 
o . 

m ■ 
O ■ 

^ : 

9 L,: 
-i— > ■ 

G ■ 
cd 

3 : 
cr 



DOI: 10. 1103/PhysRevA. 66. 062301 



I. INTRODUCTION 

Grover's search algorithm 0,0 provides an example of 
the speed-up that would be offered by quantum comput- 
ers, if and when they are built. The problem solved by 
Grover's algorithm is finding a sought-after ("marked") 
element in an unsorted database of size N. To solve this 
problem, a classical computer would need N/2 database 
queries on average, and in the worst case it would need 
N — 1 queries. Using Grover's algorithm, a quantum 
computer can find the marked state using only 0(y~N) 
quantum database queries. The importance of Grover's 
result stems from the fact that it proves the existence 
of a gap (albeit a polynomial gap) between the power 
of quantum computers and classical computers. More- 
over, the algorithm may be used to speed up the solution 
of many problems (such as NP-complete problems), for 
which no efficient classical algorithms is known. 

Along this paper we assume without loss of generality 
that N = 2", where n is an integer. The algorithm re- 
quires a register of n qubits carrying the computation. 
When we say it is in a state |x), we mean that its qubits 
are in states corresponding to the binary representation 
of the number x. Grover's original quantum search algo- 
rithm consists of the following steps: 

1. Initialize the register to H |0). That is, reset all the 
qubits to and apply the Hadamard transform to 
each of them. 

2. Repeat the following operation (named the Grover 
Iterate Q) T = times: 

(a) Rotate the marked state \k) by a phase of it 
radians 

(b) Apply the Hadamard transform to the regis- 
ter. 

(c) Rotate the |0) state by a phase of tt radians 

(d) Apply the Hadamard transform again. 

3. Measure the resulting state. 



PACS number(s): 03.67.Lx 



Several generalizations extended the original Grover al- 
gorithm. Among these is handling multiple marked 
states , and the initialization of the algorithm in any 
pure state p|. Another generalization is the replace- 
ment of the Hadamard transform by any other unitary 
operation 0, 0, 0- I n this way the algorithm may be 
used to speed-up many classical decision algorithms and 
heuristics. Other generalizations use arbitrary rotation 
angles @, replace the |0) state from Step (c) with any 
other state, or combine all of the above Q • The rotation 
angles may be tweaked in order to find a marked state 
with certainty [Tol TTll| . 

The original Grover Iterate is Q = —HI^HI^. It has 
been generalized to Q = —UI^U^IJ^, where U is an ar- 
bitrary unitary operator, s is an arbitrary state, (3 and 
7 are arbitrary angles, and M includes any number of 
marked states. We now observe that any unitary opera- 
tion Q has a unitary diagonalization. Therefore, it can 

be represented as Q = —UlgU^IJ^. This is a further 
generalization of Grover's algorithm, where the state s is 
replaced by a set of states S, each of which may have a 
different rotation angle. Thus, every iterative algorithm 
is a generalized Grover algorithm. 

In this paper we study the case where the generalized 
Grover Iterate of [fj is applied to a quantum register 
that is initialized in an arbitrary mixed state. Our study 
extends and corrects a result from ll2H. 



II. ARBITRARY PURE INITIAL STATE 

If the abovemcntioncd search algorithm is used as a 
procedure by another algorithm, it might be necessary 
to avoid its first step. Even if the initialization is per- 
formed, gate imperfection or external noise might cause 
the outcome to differ from the exact H |0) state. Rather, 
it may well be some general pure state \ipo), which is 
a superposition of the marked state and the unmarked 
states. In addition, the Iterate itslef may be imperfect: 
the Hadamard operation might be some other unitary 
operation U; the rotations of l|2a|) and (|2*c|l may be in 
angles [3 and 7 (respectively), different of 7r; and the ro- 
tated state of Step (|2*c|> may be a non-zero Is). Finally, 
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the set of sought-after items, M, may include multiple 
items. 

When the parameters of the problem are known, we 
may follow the results of Biham et al. Q, and calculate 
the probability to measure a marked state P^ as a func- 
tion of the number of Grover iterations t: 



P^o(t) = (i^o) ~ AP^ cos(2^ + 2^ ). 



(1) 



(P-ip ) and AP/, denote the average over time and the 
amplitude of P^ respectively. The subscripts denote 
that the values depend on the initial state. However, u, 
which is defined by 



E 



COS^' : : > |(i|{7|s)| COS 



+ 7 , \I-\TT\ \l 2 

—2~ +}^ \\ l \ u \ s )\ cos— ' 



is independent of the initial state. In a large search prob- 
lem with N — * 00, and with the original Grover Iterate, 
u> may be approximated by ui — -7=. With the orig- 
inal initial state H |0) studied by Grover, we re-obtain 
( p H\o)) = AP ff | ) = 5 and 0/f|o> ~ 0. 



III. ARBITRARY MIXED INITIAL STATE 

A mixed state arises when one cannot describe the 
state of a quantum system deterministically, no mat- 
ter what basis one chooses. Such a condition appears 
very often when a quantum system is entangled with 
its environment, while the environment cannot be ac- 
cessed or manipulated. The state of such a system may 
be described by a completely-positive trace- 1 hermitian 
density matrix, denoted by p. An equivalent descrip- 
tion is an ensemble £ — {p^, where = 1 
and p — X^fV IVv) (VVil- According to this description, 
the system is in the pure state \ipn) with probability p^. 
When a unitary operation V is applied to the mixed state 
it transforms the state into {%\V^. The 
mixedness of the state does not change, and it may be 
thought as if V transforms each of the components of £ 
independently of the others. 

Extending the argument of Section [HI the initial state 
of the quantum register might not be pure, due to ex- 
ternal noise, decoherence or previous manipulations. In- 
stead, the initial state may be some general mixed state 
£. Given the description of £ as an ensemble, all we 
can say is that the register is in the pure state with 
probability p^ (for all i's). 

When the Grover algorithm is applied to a register 
which is in a pure state |Vv)> the probability to measure 
the marked state is P M (t) . The probability for the register 
to be in \tp^ is p^ (Considering the ensemble £). Thus, 
the total probability to measure the marked state is the 
weighted average 



P(t) = 5> M P M (t) 
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FIG. 1: The change of probabilities as projections of rotating 
vectors. 



P, ((^ 



AP M cos 



(2^ + 2^)). (2) 



The functions 



P»{t) - (P„) = -AP M cos {2ut + 20„) 

share a sinusoidal form, differing in amplitude and phase, 
but not in frequency. They may be thought of as the 
projections of vectors rotating in frequency u>, as exem- 
plified by Figure ^ Therefore, their weighted sum AP 
(the Center of Mass of the vectors in the figure) is a si- 
nusoidal function with the same frequency: 



P(t) = (P) - AP cos ( 2ut + 2<j) 



(3) 



while 



(P)=5> M <P M ), 



AP 



and 



\ 



^2 P^P^ cos 2(j)A + E P^ P n sin 



tan 2(f) 



E^AP^sin^) 

Eu-Pm ap m cos ( 2( /v) 



(4) 



(5) 



The probability to measure a marked state reaches its 
maximum value 



after T 



2uj 



P max = (P) + AP 



iterations. 
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If the algorithm is repeated until success with T iter- 
ations each time, the expected total time to measure a 
marked state is 



Tn = 



7T - 26 



Pmax 

since the number of repetition until success is distributed 
geometrically with parameter P ma x- When the original 
Iterate is used, and a single item is sought after, this 

reduces to Td = ™~ \fN. If this value is significantly 
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smaller than the classical expected time Tq< = N/2, then 
the quantum algorithm has an advantage. Quantita- 
tively, the expected number of oracle queries that the 
quantum algorithm requires is smaller by a factor of 



Tc 
To 



NljP„ 



2(f> ir -2<f> 



IV. EXAMPLES 



(6) 



For clarity and simplicity, our examples use the original 
Grover Iterate and single marked state | k) , with different 
initial mixed states. 



Pure Initial State 



C. Initial State Where m of the Qubits Are Mixed 

Let us study the case where the register is initialized 
to p m -mix = Y^i^ 1 H \i) (i\ H. This state may oc- 
cur if the m least significant qubits of the register are 
totally mixed before the first Hadamard transform is ap- 
plied. Since all H\i) are orthogonal to P|0) (except for 
H |0) itself) and they are almost orthogonal to |fc) (since 
|(fc| H |i)| 2 = jr), the evolution of p m -mix is governed by 
{p = 2- m , H |0)} and we obtain (P) = AP = 5^ and 
6 = 0. Large m would render the algorithm useless. 



V. ALGORITHM USEFULNESS AND 
ENTROPY 

The von Neumann entropy of a mixed state p is de- 
fined as S(p) = — trplogp. Bose et al. presented a 
new model for quantum computation and laid out a new 
proof for the optimality of the Grover algorithm. How- 
ever, one of their results was the following: if the Grover 
algorithm is initiated with a mixed state p, such that 
S(p) > \ log N, the algorithm would have no advantage 
comparing to the classical case. This is in disagreement 
with our findings A counter-example to their claim 
is p 1 as defined above. The entropy of pseudo- 

log N " 

pure state is 



When the arbitrary mixed state is chosen to be pure, 
the summations are degenerated and the results of |4| are 
achieved. For example, if the initial state is the original 
£ = {p = l,H\0)}, the original Grover case is found. 

If 8 = {p = 1, |fc)}, then (P) = AP = \ and 4> = f . 
An interesting known property of the Grover algorithm 
is that for all states orthogonal to both \k) and H |0), 
(P) = AP = 0. 



B. Pseudo-Pure Initial State 

Ensembles where a state \4>) appears with probability 
e + and any state orthogonal to it appears with equal 
probabilities of are called pseudo-pure mixed states. 
They are written more conveniently as p e - P ure = (1 — 
e)^ + e (tp\. Notice that < e < 1 is a measure of 
the purity of p: when e = it is totally mixed, and when 
e = 1 it is totally pure. It is easy to see that in the limit 
of large TV, (P) = e (P v ,), AP = eAP^ and 4> = For 
example, for 



1 



)4 + 



1 



logN'N \ogN 



H |0> (0| H, 



we obtain (P) = AP = 



and 6 — 0. Notice that 



21ogAf 

although p is extremely mixed, the quantum advantage 



is of factor 



2y^V 

TT log N ' 



$ \Pe-pure) 



AT-1 



E V log 



e|0) <0| 
1 - e 



N 



l + (JV-l)c, l + (JV-l)e 
-log- 



N 

-(TV-i)i-^log 1 ' 



N 



N 



N 



l + (N-l)e, l + (N-l)e 
log , 

N N 

and for large N, where Nj (N — 1) w 1, 

- -(l-e)]og^-(l +e ) 108(1 
= (l-e)logTV-(l-e)log(l-e) 

^ + £ ) l0g (^ 
(l-e)logAT-f, 



(7) 



where I = (1 - e) log(l - e) + + e) log (j? + e) € 
(-1, 0.8) for any < e < 1 and any N > 2. For e = 



we obtain 5(p_i__ p „ re ) = [l - logiV + 0(1) = 

log N + 0(1). This entropy is almost maximal. However, 
as noted above, the Grover algorithm outperforms any 
classical algorithm, even when it is initialized with this 
state. 
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Entropy is not a good measure for the usefulness of 
Grover's algorithm. For practically every value of en- 
tropythere exist states that are good initializers and 
states that are not. For example, S(p( n _iy mix ) — 
logiV — 1 = S(p ^ i pu rc )- but when initialized in 
P(n-i)-mix; the Grover algorithm is as bad as guessing 
the marked state. Another example may be given us- 
ing the pure states H |0) (0| H and H\l) (1| H. With the 
first, Grover arrives to the marked state with quadratic 
speed-up, while the second state is practically unchanged 
by the algorithm. 



VI. ACKNOWLEDGMENTS 



We thank Tal Mor for valuable discussions and sug- 
gestions that made the compiling of this paper possible. 
The work was partially supported by the European Com- 
mission through the 1ST Programme under contract IST- 
1999-11234. The first author was partially supported by 
the fund for the promotion of research at the Technion 
and the Israel MoD Research and Technology Unit. 



[1] L. K. Grover, in Proceedings of the Twenty-Eighth Annual 

ACM Symposium on Theory of Computing (ACM Press, 

New York, 1996), pp. 212-219. 
[2] L. K. Graver. Phys. Rev. Lett. 79. 325 (1997), also hi 

quant-p h/9706033 
[3] M. Boyer, G. Brassard, P. H0yer, and A. Tapp, Fortsch. 

Phys. 46, 493 (1998), also in quant-ph/9605034 
[4] E. Biham et ai, Phys. Rev. A 60, 2742 (1999), also in 

quant-ph/9711070 
[5] L. K. Grover, Phys. Rev. Lett. 80. 4329 (1998), also in 

quant-ph/9712011 
[6] R. M. Gingrich, C. P. Williams, and N. J. Cerf, Ph ys. 

Rev. A 61, 052313 (2000), also in quant-ph/9904049 
[7] G. Brassard, P. H0yer, M. Mosca, and A. Tapp, 

Quantum Amplitude Amplification and Estimation, In 

quant-ph/0005055 2000. 
[8] G. L. Long, Y. S. Li, W. L. Zhang, and C. C. Tu, Phys. 



Rev. A 61, 042305 (2000). 
[9] E. Biham et ai, Phys. Rev. A 63, 012310 (2001), also in 
|quant-ph/0010077| 

[10] P. H0yer, Phys. Rev. A 62, 052304 (2000). 

[11] G. L. Long, L. Xiao, and Y. Sun, General Phase Matching 
Condition for Quantum Searching, In quant-ph/0107013 

[12] S. Bose, L. Rallan, and V. Vedral, Phy s. Rev. Lett. 85, 
5448 (2000), also in quant-ph/0003072 

[13] Bose et al. have a mistake regarding the entropy of the 
classical search problem. Just above their Equation (10), 
they say that classical search can change entropy by 
log \/N in \/ r N steps. This is true for a search field of 
size VN , but wrong for the question in matter where the 
search field is of size N . This mistake does not invalidate 
the other results in their paper. 



